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1  Introduction 


Computer  programs  are  now  used  routinely  to  support  a  wide  variety  of  ship-related 
research.  At  DREA,  computer  programs  have  been  developed  to  aid  in  ship  design  as  well  as 
to  determine  ship  motions,  ship  structural  vibrations,  and  flows  around  the  hull  and  into  the 
propeller  plane.  Each  of  these  programs  requires,  as  input,  a  description  of  the  geometry  of 
the  ship  hull,  but  each,  at  present,  requires  a  different  format  for  this  description.  Thus,  when 
performing  a  series  of  calculations  on  the  same  hull,  much  time  is  wasted  preparing  several 
different  data  files  each  of  which  describes  the  hull  geometry.  To  circumvent  this  problem,  a 
system  for  describing  the  hull  geometry  has  been  developed  which  is  capable  of  describing  all 
relevant  characteristics  of  the  hull  for  the  above  programs.  The  system  is  called  HLLSRF. 

The  HLLSRF  hull  representation  system  has  grown  out  of  earlier  hull  modelling  programs 
associated  with  the  HLLFLO  programs  for  calculating  the  flow  around  hulls1,2,3,4.  HLLSRF 
has  extended  these  programs  to  allow  the  representation  of  many  common  hull  features  which 
were  not  allowed  previously:  notably,  bulbous  bows,  complex  stern  profiles,  slanted  transoms, 
and  stepped  decks. 

A  useful  feature  of  the  HLLSRF  representation  scheme  is  that  variations  of  existing  hulls 
can  be  represented  by  specifying  how  they  differ  from  their  parent.  For  example,  the  length 
of  the  hull  could  be  changed  simply  by  specifying  the  new  length  in  a  modification  file;  all 
other  aspects  of  the  hull  would  remain  unchanged.  Similarly,  the  position  of  the  propeller  disk 
could  be  altered  or  its  draft  lowered.  Using  a  method  similar  to  that  of  Lackenby5,  the  form 
of  the  hull  can  also  be  modified  to  change  its  prismatic  coefficients  or  its  longitudinal  centre  of 
buoyancy. 

The  main  function  of  a  HLLSRF  representation  is  to  provide  a  regular  coordinate  system 
covering  the  entire  hull.  The  coordinate  system  provides  the  means  for  performing  a  wide 
variety  of  hydrodynamic  calculations  ranging  in  complexity  from  simple  evaluation  of  sectional 
properties  for  strip  theory,  to  the  solution  of  complicated  differential  equations  to  determine 
the  viscous  flow  around  the  hull.  HLLSRF  has  been  designed  to  maintain  the  regularity  of  this 
coordinate  system  so  that  numerical  difficulties  can  be  avoided.  To  facilitate  flow  calculations, 
an  additional  regular  coordinate  system  is  provided  which  covers  only  the  portion  of  the  hull 
below  the  waterline.  Using  the  method  of  hull  modification  described  above,  the  waterline  itself 
may  be  changed  without  changing  the  rest  of  the  hull  representation.  This  allows  a  simple  and 
efficient  means  for  calculating  flows  at  different  drafts,  at  different  trim  angles,  or  at  different 
Froude  numbers. 

HLLSRF  hull  representations  are  also  extensible:  extra  data  may  be  included  in  the  data 
file  describing  a  hull.  To  make  use  of  the  extra  data,  the  user  must  supply  subroutines  which 
read  and  access  the  data  in  a  manner  appropriate  to  his  needs;  otherwise  they  will  be  ignored. 
The  extensibility  of  HLLSRF  is  very  convenient,  since  it  allows  for  the  inclusion  of  data  which 
may  be  required  for  one  program  but  not  for  others.  It  also  allows  a  convenient  means  for 
extending  the  system  to  include  hull  features  not  considered  important  now  but  which  may  be 
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2  Hull  Representation  Strategy 


In  this  section  an  overview  of  the  HLLSRF  representation  is  given.  A  more  precise  defini¬ 
tion  is  given  in  Section  3. 


2.1  Coordinate  Systems 

The  hull  is  defined  in  terms  of  a  Cartesian  coordinate  system  whose  axes  are  aligned  with 
the  waterlines  and  centreplane.  The  three  Cartesian  coordinates  are  denoted  (X,Y,Z)  and  are 
defined  as  follows. 

•  The  coordinate  X  increases  from  bow  to  stern.  The  forward  perpendicular  is  X  —  0  and 
the  aft  perpendicular  is  X  =  L,  where  L  is  the  length  between  perpendiculars.  The  planes 
of  constant  X  are  hull  sections. 

•  The  coordinate  Y  increases  from  the  centreplane  to  the  point  of  maximum  breadth.  The 
centreplane  is  Y  =  0.  Owing  to  the  symmetry  of  the  hull,  it  is  usually  not  necessary 
to  decide  whether  Y  increases  to  port  or  to  starboard,  but,  to  maintain  a  right-handed 
coordinate  system,  Y  should  increase  to  starboard.  Surfaces  of  constant  Y  are  parallel 
to  the  centreplane. 

•  The  coordinate  Z  increases  from  the  baseline  to  the  point  of  maximum  height  above  the 
baseline.  The  surfaces  of  constant  Z  are  waterplanes.  The  waterplane  which  contains  the 
baseline  is  Z  —  0. 

The  position  vector  (X,Y,Z)  will  sometimes  be  denoted  by  the  single  bold  face  symbol  X. 

A  hull  representation  is  simply  a  parameterization  of  the  hull  surface  so  that  the  Cartesian 
coordinate  vector,  X,  of  a  point  on  the  hull  can  be  obtained  as  functions  of  the  two  parameters 
(x,y).  That  is, 

X  =  f(x,y)  (2.1) 

for  some  vector  function  f  =  (fx,  fy >  fz)-  The  parameters  x  and  y  define  a  non-Cartesian 
coordinate  system  on  the  hull  surface.  In  general  x  increases  along  the  length  of  the  hull  (lines 
of  constant  y  run  from  bow  to  stern)  while  y  increases  around  its  girth  (lines  of  constant  x  run 
from  the  keel  to  the  deck  edge).  Both  coordinate  systems  are  summarized  in  Figure  1. 

For  computer  applications,  the  vector  function  f  is  usually  expressed  as  a  linear  combina¬ 
tion  of  a  series  of  basis  functions,  bn(x )  and  bj(y). 

N,  Ny 

«jn&j {y)bn(x)  (2.2) 

n=l j  =  1 

Each  hull  is  then  described  by  the  coefficients  an  which  may  be  stored  in  a  data  file. 
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Figure  2:  A  ship  hull  with  a  stepped  deck 


Figure  3:  xy-coordinate  grid  which  preserves  a  rectangular  domain 


Forward  x — »  Aft 


Figure  4:  Non-rectangular  xy-domain  used  to  generate  Figure  2 


\ _ \ 

Figure  6:  Surface  S  for  Figure  5 


Figure  7:  (x,y)  domain  for  Figure  5 
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Figure  8:  Hull  with  a  chine 


Figure  9:  Stern  with  a  propeller  bossing 
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which,  in  the  way  it  is  represented,  is  essentially  the  same  as  the  chine  line  of  Figure  8;  the 
function  df/dy  is  discontinuous.  The  second  is  the  corner  where  the  stern  frame  meets  the  flat 
of  keel;  in  this  case  the  function  df/dx  is  discontinuous.  The  discontinuity  is  localized  to  the 
region  of  the  keel. 
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otherwise  unconstrained.  When  it  is  important  to  be  able  to  model  the  hull  on  both  sides  of 
the  centreplane,  y  is  allowed  to  vary  between  —1.0  and  1.0. 

X(x,y)=  fx(x,-y)  ' 

Y(x,y)=-fy(x,-y)  -  if  —  1  <  y  <  0  (3.2) 

Z(x,y)=  fz(x,-y )  _ 

Before  the  functional  forms  of  fy  and  fz  can  be  specified,  it  is  necessary  to  choose  a  set 
of  basis  functions  bn(x,y).  HLLSRF  uses  B-splines  of  arbitrary  order  (for  further  information 
on  B-splines  see  de  Boor6).  B-splines  are  a  popular  choice  as  the  basis  functions  for  hull 
representations2,7-16  as  they  can  be  calculated  both  accurately  and  efficiently  and  are  capable 
of  describing  a  wide  variety  of  complex  shapes7'8,9,10.  By  manipulation  of  the  B-spline  knots, 
it  is  possible  to  generate  the  ‘cuts’  in  the  hull  and  the  discontinuities  in  derivative  described  in 
Section  2.4;  details  are  given  in  Section  3.3. 

HLLSRF  uses  the  following  functional  forms  to  define  the  Cartesian  coordinates  Y  and  Z 
in  terms  of  B-splines. 


N»  Nx 


Y  =  fy(x,y)  =  ByyajnBnikx(x)BJtky(y) 

(3.3) 

j= ln=l 

Ny  Nx 

Z  =  fz(x,y)  =  D'£,'£fijnBnM,{x)Bjjk,(y) 

(3.4) 

j=l n=l 


where 

B  is  the  half-breadth  of  the  ship, 

D  is  its  depth  from  the  baseline  to  the  highest  point  above  the  baseline, 

Bn<kx(x)  is  a  B-spline  of  order  kx  corresponding  to  a  knot  sequence  t\x\  i  =  1, . . . ,  Nx  +  kx,  for 
each  n  =  1, . . . ,  Nx, 

Bn,kv(y )  is  a  B-spline  of  order  ky  corresponding  to  a  knot  sequence  t\v\i  =  1, . . . ,  Ny  +  ky,  for 
each  n  =  1,. . . , Ny,  and 

are  coefficients, 

The  B-splines  are  normalized  to  sum  to  one: 

N 

T  Bi,k{x)  =  1  if  tk  <  x  <  tu-k+1  (3.5) 

«=i 


y  € 


The  sums  of  the  B-splincs  return  non- zero  values  for  Y  and  Z  if  x  6  (<ir\  t ^l+kx )  an(i 
(*!»*%+*„)•  Hence  it  is  necessary  that  t[x)  <  xmtn,  t^]+kx  >  xmax ,  t[v)  <  0  and 


>  1. 
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derivative  of  yd(x );  it  may  be  generated  by  including  a  knot  of  multiplicity  k  —  1.  The  number 
of  deck  edge  discontinuities,  the  station  numbers  at  which  they  occur,  and  their  types  (i.e. 
whether  they  are  steps  or  corners)  are  all  stored  in  additional  HLLSRF  data  structures. 

3.1.2  Bow  and  Stern  Profile  Curves 

The  functions  y™(x)  used  to  approximate  the  bow  profile  curve  are  also  specified  in  terms 
of  B-splines.  Thus,  for  given  x,  the  point  ( x,y )  lies  (approximately)  on  the  bow  profile  curve  if 
am  <  x  <  bm  (see  equation  (2.3))  and  y  =  y™(x)  for  some  m  €  [1,M{,]  where 

»”(*)  =  £«■%-(*)  <3-8> 

n=  1 

and 

B™km(x)  is  a  B-spline  of  order  km  corresponding  to  a  knot  sequence  tj",  i  =  1, . . . ,  Nm  +  km , 
for  each  n  =  1, . . . ,  Nm,  and 

f 3 ™  are  coefficients, 

At  the  stern  a  completely  analogous  set  of  splines  is  used  to  represent  the  stern  profile 
curve. 


3.2  Slanted  Transoms 

HLLSRF  represents  a  transom  as  a  plane  which  intersects  the  hull  near  the  stern  (see 
Figure  12).  It  does  not  allow  the  representation  of  curved  transoms  or  more  complex  transoms 
consisting  of  several  different  planes  as  in  Figure  13;  however,  for  most  purposes,  the  hull  of 
Figure  13  could  be  adequately  represented  by  including  the  line  A-B  in  the  deck  edge  rather 
than  the  transom. 

The  transom  plane  is  defined  by  specifying  a  point  through  which  it  passes,  and  an  aft- 
facing  normal  vector.  Let  the  Cartesian  coordinate  vector  of  the  point  through  which  the 
transom  passes  be  X*  and  its  normal  n.  A  point  X  lies  on  the  transom  if 

(X  —  X()  •  n  =  0  (3.9) 

By  symmetry,  ny  is  zero,  and  one  can  choose  X<  to  lie  on  the  centreplane  so  that  \\  =  0. 
Hence,  to  define  the  transom,  only  Xt ,  Zt,  nx  and  nz  need  be  given. 

When  a  transom  is  defined,  the  hull  domain  V  must  be  restricted  to  discard  points  which 
lie  beyond  the  transom:  V  =  {(ar,y)  :  (X(x,y)  -  X*)  •  n  <  0}. 


3.3  Hull  Discontinuities 

Discontinuities  may  be  accommodated  in  B-spline  hull  representations  in  three  different 
ways.  Rogers  and  Satterfield11,  Kukner  and  Ergur12  and  Rosovic13  generate  knuckles  by  a 
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simple  manipulation  of  the  hull  B-spline  coefficients  ajn  and  /3jn  (see  equations  (3.3)  and 
(3.4)).  If  Oj— fc-f iftl  —  Otj—k+ 2,n  =  ...  =  Qj—  l,n  3.nd  /3j—k+l,n  —  0j  —  k+ 2,n  ~  =  0j—l,n  for  «H1 
n  €  [nm,n,nmar],  then  the  line  y  =  tj  is  a  knuckle  for  £  €  [^nm.ri+fcI,<nmaI].  However,  Hally14 
has  pointed  out  that  the  xy-coordinate  generated  by  this  method  is  singular  at  the  knuckle; 
the  singularity  will  often  cause  difficulties  when  used  in  numerical  calculations.  Moreover,  in 
HLLSRF,  because  the  Cartesian  coordinate  X  is  not  described  using  B-splines,  this  method 
cannot  be  used  to  generate  hull  discontinuities  along  a  hull  station  (as  in  Figure  10,  for  example). 
Furthermore,  this  method  cannot  be  used  to  generate  complete  discontinuities  in  the  surface  S 
such  as  the  ‘cut’  used  to  define  the  propeller  bossing  in  Figure  10.  Hence,  a  better  method  for 
representing  hull  discontinuities  is  required. 

The  other  two  methods  for  generating  knuckles  use  the  fact  that  multiple  knots  can  be 
used  to  decrease  the  degree  of  continuity  of  a  spline  (see  de  Boor6,  chapter  IX);  in  particular, 
if  a  knot  sequence  of  a  B-spline  of  order  k  has  a  knot  of  multiplicity  u,  then  the  spline  curve 
can  have  discontinuous  ( k  —  i/)th  derivative  at  that  knot.  Knuckles  can  be  defined  by  having 
knots  of  multiplicity  ky  —  1  in  the  knot  sequence  t\vK 

Hayes10,  Stroobant  and  Mars9,  and  Hally3  have  generated  knuckles  by  allowing  the  knot 
sequence  to  vary  with  x.  Suppose  that  { t\y is  such  that  all  knots  are  distinct  for 

x  <  xio  and  x  >  £/,,,  but  that  A:y  —  1  knots  coalesce  when  x/0  <  x  <  x^,-.  Then  the  hull  will 
generally  exhibit  a  discontinuity  in  derivatives  with  respect  to  y  for  xi0  <  x  <  xm  but  will  be 
smooth  outside  this  range;  that  is,  there  will  be  a  knuckle  beginning  at  x/0  and  extending  to 
Xhi ■  Similarly,  discontinuities  along  hull  stations  can  be  defined  by  letting  the  knot  sequence 
vary  with  y. 

While  this  method  seems  to  work  well  when  defining  the  hull  representation,  it  has  dis¬ 
tinct  shortcomings  when  calculations  are  performed  using  the  representation.  In  particular, 
evaluation  of  points  on  the  hull  given  ( x,y )  is  less  efficient  since  the  knot  sequence  must  be 
calculated  before  the  B-splines  can  be  evaluated  and  summed.  Techniques  for  speeding  up  the 
calculation  are  described  by  Hally3.  The  inefficiency  is  much  worse  when  y- derivatives  of  the 
functions  Y(x,y)  and  Z(x,y)  are  to  be  calculated  (for  example,  when  calculating  a  normal  to 
the  hull  at  a  given  point)  since  now  the  derivatives  of  the  B-splines  with  respect  to  the  knots 
must  also  be  calculated.  This  overhead  is  too  costly  for  hull  representations  which  are  to  be 
used  in  hydrodynamic  computations;  it  is  not  used  by  HLLSRF. 

The  third  method  also  uses  multiple  knots  to  generate  the  knuckle  but  does  so  while 
keeping  the  knot  sequences  fixed.  A  knot  of  multiplicity  v  in  the  knot  sequence  does  not 
necessarily  imply  that  the  spline  will  have  a  discontinuous  ( k  —  i/)th  derivative  at  the  knot; 
it  merely  allows  the  spline  to  a  have  a  discontinuous  ( k  -  i/)th  derivative.  In  fact,  by  proper 
adjustment  of  the  spline  coefficients,  the  spline  can  be  kept  smooth  at  the  multiple  knot.  The 
third  method  of  knuckle  representation  generates  a  discontinuity  in  the  spline  where  the  knuckle 
is  intended,  but  adjusts  the  spline  coefficients  so  that  the  splines  are  smooth  elsewhere. 

For  example,  suppose  t\y^  =  ...  =  =  f,  so  that  there  is  a  knot  of  multiplicity  ky  —  1 

at  y  =  t.  The  hull  will  then  normally  have  discontinuous  first  derivative  with  respect  to  y 
at  y  =  t:  i.e.  there  will  be  a  knuckle.  By  adjusting  the  spline  coefficients  ajn  and  (5]n  for 
n  =  r, . . .  ,s,  one  can  remove  the  discontinuity  over  the  range  x  €  ^l*x+j].  Hence,  by 
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subsidiary  spline  curve,  y  =  yw(x). 


Vw(x) 


N 

^2  PnBn,k(x)  for  a<  x  <b 

n=  1 

0  otherwise 


(3.10) 


If  the  waterline  curve  is  not  present  in  the  hull  representation,  it  is  assumed  that  the  waterline 
curve  is  y  =  1.  Use  of  the  waterline  curve  avoids  regeneration  of  the  hull  whenever  the  waterline 
is  changed:  for  example,  if  there  is  a  change  of  draft  or  Froude  number.  However,  the  waterline 
will  not  usually  be  a  line  of  constant  y ,  a  property  which  is  desirable  for  performing  calculations 
below  the  waterline.  Therefore  a  new  coordinate,  s,  is  defined  by 


s  = 


y 

Vw(x) 


(3.11) 


It  is  0.0  at  the  keel  and  1.0  at  the  waterline.  The  coordinate  system  defined  by  (x,s)  is  preferred 
for  most  hydrodynamic  calculations.  The  Cartesian  coordinates  corresponding  to  a  point  (x,s) 
are  easily  obtained  using  the  inverse  relation  y  =  syw(x). 

X(x,s)  =  t(x,syw(x))  (3.12) 


Derivatives  may  be  obtained  using  the  chain  rule.  Use  of  (x,s)  coordinates  implies  some  loss 
of  efficiency  in  comparison  with  (x,y),  but  this  is  easily  regained  by  the  requirement  that  only 
one  hull  need  be  generated  to  accommodate  all  possible  waterlines.  The  time  taken  to  calculate 
the  spline  representation  of  the  waterline  curve  yw(x)  is  not  significant. 


3.4.2  The  Curve  of  Sectional  Areas 


Many  strip  theory  programs  use  sectional  properties  of  the  hull  that  can  be  determined 
once  the  curve  of  sectional  areas  is  known.  To  avoid  repeated  calculation  of  the  sectional 
areas,  the  HLLSRF  representation  allows  the  representation  of  the  curve  of  sectional  areas  by 
a  subsidiary  spline. 

(  N 


A{x) 


^2  f3nBn,k(x )  for  o  <  x  <  b 

n=l 


(3.13) 


0 


otherwise 


where 


A(x)  is  the  area  of  the  half-section  at  x  divided  by  BZdwi 
B  is  the  half-breadth  of  the  ship,  and 


Zdw  is  the  design  waterline  of  the  ship. 
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4  Hull  Modifications 


One  of  the  most  useful  features  of  the  HLLSRF  representation  scheme  is  that  variations 
of  existing  hulls  can  be  generated  using  modification  files:  that  is,  a  new  hull  can  be  defined 
by  specifying  how  it  differs  from  an  existing  hull.  The  changes  are  specified  by  data  in  the 
modification  file.  For  example,  the  waterline  for  a  hull  may  be  changed  by  specifying  a  new 
spline  representation  for  the  waterline  curve  in  a  modification  file.  All  other  aspects  of  the  hull 
remain  unchanged.  Similarly,  the  position  of  the  propeller  disk  could  be  altered  or  the  deck 
edge  changed. 

Changes  to  the  hull  itself  are  also  possible.  The  simplest  of  these  involve  changes  of  the 
hull  dimensions,  L ,  B,  and  D.  The  form  of  equations  (3.1),  (3.3)  and  (3.4)  has  been  chosen  so 
that  these  quantities  can  be  changed  without  the  need  of  altering  the  spline  coefficients. 

Lackenby5  first  introduced  the  idea  of  modifying  hulls  by  shifting  their  stations  between 
the  perpendiculars.  This  idea  has  been  developed  for  use  in  computer  hull  representations  by 
Hally16  and  Hally  and  Gilbert4.  The  shift  of  stations  is  specified  by  defining  a  transformation 
between  the  station  numbers  on  the  modified  hull  and  the  station  numbers  on  the  parent.  The 
modified  hull  is  defined  by 


(z  -  xppm)L 

(4.1) 

xAPm  ~  xFPm 

fy  (xp(x)) 

(4.2) 

fz(xP(x)) 

(4.3) 

where 

xppm  is  the  station  number  of  the  forward  perpendicular  of  the  modified  hull, 

XAPm  is  the  station  number  of  the  aft  perpendicular  of  the  modified  hull,  and 

xp(x)  is  a  function  which  returns  the  station  number  on  the  parent  hull  which  corresponds  to 
the  station  x  on  the  modified  hull. 

When  the  hull  is  modified  in  this  way,  there  is  no  need  to  modify  any  of  the  subsidiary  splines 
in  the  hull  representation.  For  example,  the  waterline  of  a  hull  with  modified  station  numbers 
is  generated  from  that  of  the  parent  hull  by 

V  =  yw{Xp(x))  (4.4) 

The  function  xp(x)  is  represented  by  B-splines  as  follows. 

N 

xp(x)  =  xp(x)+  Y^0nBnAxP(x))  (4-5) 

n=l 
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5  Using  the  HLLSRF  Library 


The  HLLSRF  library  consists  of  many  sub-programs  which  perform  calculations  using 
the  HLLSRF  hull  representation  of  a  hull.  The  most  commonly  used  calculations  for  which 
sub-programs  are  provided  are  described  in  this  section. 


5.1  Loading  the  Representation  of  a  Hull 


Data  defining  a  HLLSRF  hull  representation  is  stored  in  a  computer  file  organized  into  a 
series  of  labelled  records  using  the  format  described  by  Hally17;  details  of  the  records  recognized 
by  HLLSRF  are  described  in  the  HLLSRF  User’s  Guide18. 

Before  any  calculations  can  be  performed,  a  hull  representation  must  be  loaded  by  reading 
a  data  file.  Although  there  are  two  two  different  types  of  hull  data  files,  one  for  a  parent  hull 
and  one  for  hull  modifications  (see  Section  4),  a  single  sub-program  is  used  to  read  all  files:  the 
existence  of  modifications  is  transparent  to  the  user.  This  is  also  true  for  all  of  the  HLLSRF 
sub-programs  described  here:  they  all  work  equally  well  whether  the  hull  has  been  modified  or 
not. 


5.2  Coordinate  Transformations 

As  was  pointed  out  in  the  introduction,  the  main  purpose  of  the  hull  representation  is  to 
provide  coordinate  transformations  between  the  Cartesian  coordinates  (X,Y,Z)  and  the  non- 
Cartesian  hull  coordinates  (x,j/).  Accordingly,  HLLSRF  provides  sub-programs  for  performing 
such  coordinate  transformations. 

A  sub-program  is  provided  which,  given  a  point  (x,y),  will  calculate  any  of  the  Cartesian 
components  X,  Y ,  and  Z,  or  any  of  their  derivatives  with  respect  to  x  and  y.  This  involves  a 
simple  evaluation  of  the  splines  according  to  equations  (3.1)  (3.3)  and  (3.4).  When  a  modifi¬ 
cation  of  hull  form  has  been  specified  (see  Section  4),  then  equations  (4.1),  (4.2),  (4.3),  (4.6), 
and  (4.5)  are  used. 

For  hydrodynamic  calculations,  it  is  more  convenient  to  use  the  coordinates  (x,s)  which 
are  defined  only  on  the  submerged  portion  of  the  hull  (see  Section  3.4.1).  A  sub-program  is 
provided  which,  given  a  point  (x,y),  will  calculate  any  of  the  Cartesian  components  X ,  Y,  and 
Z,  or  any  of  their  derivatives  (up  to  second  order)  with  respect  to  x  and  s.  This  involves  simple 
evaluation  of  the  hull  splines,  the  waterline  spline,  ^(x),  and  their  derivatives. 

The  HLLSRF  library  also  provides  sub-programs  which  will  determine  points  (x,y)  or  (x,s) 
given  values  of  the  Cartesian  components.  Given  a  value  for  x  and  a  value  for  Y,  say  Yo,  there 
is  a  sub-program  which  will  determine  the  value  of  y  (or  s)  such  that  fy{x,y)  =  Y0-  Similarly, 
a  value  of  Z  could  be  specified.  Likewise,  if  y  (or  s)  is  specified,  there  is  a  sub-program  which 
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a  vector  to  which  it  is  tangent,  t.  The  line  can  then  be  parameterized  in  terms  of  the  variable 
u  by:  {X  :  X  =  Xo  +  ut}.  A  sub-program  is  provided  which  will  determine  a  point  (x,s)  at 
which  the  line  intersects  the  hull:  that  is,  it  solves  the  three  simultaneous  equations 

X0  +  tit  =  X(x,y)  (5.1) 

for  the  unknowns  x,  y,  and  u.  For  example,  by  setting  Xq  to  the  centre  of  the  propeller  disk 
and  t  to  its  normal,  this  sub-program  could  be  used  to  determine  the  point  at  which  a  propeller 
shaft  meets  the  hull. 

A  plane  may  be  defined  by  specifying  a  Cartesian  point  through  which  it  passes,  Xo,  and 
its  normal,  n.  A  point  ( x,y )  lies  on  the  intersection  of  this  plane  and  the  hull  if 

(X(x,y)  -  X0)  •  n  =  0  (5.2) 

Sub-programs  are  provided  which  will  solve  this  equation  for  x  given  a  value  for  y,  or  which  will 
solve  it  for  y  given  a  value  for  x.  In  addition,  a  sub-program  is  provided  which  will  calculate 
a  series  of  points  on  the  line  of  intersection  of  a  plane  with  the  hull.  This  could  be  used,  for 
example,  to  determine  the  intersection  of  the  waterplane  with  the  hull  for  different  trim  or  roll 
angles. 
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1.  The  values  of  y"n  are  set  to  the  fractional  arclength  around  the  station,  where  the  ar- 
clength  is  approximated  by  straight  line  segments  joining  the  offset  points. 


Vln 


=  0 


Vmn 


E£2  y/&n  -  Yj- l,n)2  +  (Zin  -  Z,-i,n)2 
E&  \J(Y>n  ~  Yi-  1>n)2  +  (Zin  ~  ^,-l.n)2 


(6.1) 

(6.2) 


This  method  was  originally  proposed  by  Sachdeva  and  Preston20  and  was  used  as  the 
basis  of  the  hull  representation  system  in  the  DREA  hydrodynamic  package  HLLFLO1,2. 
The  fractional  arclength  method  has  the  advantage  that  y"  n  is  always  well-defined  for  any 
set  of  non-coincident  offset  points;  however,  it  has  the  disadvantage  that  the  value  of  y!^n 
depends  on  all  the  offset  points  on  the  station.  Adverse  effects  from  this  are  illustrated 
in  Figure  14  which  shows  the  y"-coordinate  lines  on  the  afterbody  of  a  hull  with  a  skeg 
(and  no  knuckles).  The  corner  in  the  skeg  at  the  keel  causes  a  ‘kink’  in  all  the  lines  of 
constant  y"  all  the  way  to  the  deck  edge,  although  with  severity  which  decreases  as  the 
deck  edge  is  approached.  Thus  local  variations  in  geometry  can  have  non-local  effects  in 
the  definition  of  y" . 

2.  The  y^n  are  set  to  the  fractional  angle  subtended  by  the  point  at  the  intersection  of  the 
centreplane  and  the  deck.  This  method  requires  that  Kin  =  0. 


Vmn  =  -arctan 


2 
— i 
7T 


(  Knn  \ 

\  zMn  n  —  Zmn  J 


(6.3) 


The  fractional  angle  method  has  been  used  by  Muraoka  and  Shirose  in  calculations  of 
turbulent  boundary  layer  flow  around  ship  hulls21. 

If  ZMn<n  is  the  same  on  all  stations,  as  is  often  the  case,  the  fractional  angle  method  has 
the  advantage  that  y"„  depends  only  on  the  offset  point  (Ymn,Zm „);  thus,  the  non-local 
effects  found  with  the  fractional  arclength  method  can  be  avoided.  Figure  15  shows  the 
y"-lines  generated  for  the  same  hull  as  in  Figure  14. 

On  the  other  hand,  this  method  does  not  work  well  on  stations  which  are  long  and  thin,  as 
often  occur  near  the  bow.  Figure  16  shows  the  distribution  of  y"-values  on  such  a  station 
using  this  method.  The  clustering  of  the  y"-lines  at  the  upper  part  of  the  station  is  not 
desirable.  Moreover,  in  the  limit  that  the  width  of  the  station  vanishes,  this  method  of 
calculation  fails.  It  will  also  fail  when  there  are  offset  points  with  negative  Y-values,  for 
then  it  would  be  possible  to  have  two  different  offset  points  with  the  same  value  of  y". 

3.  The  cartesian  components  Y  and  Z  are  first  normalized  by  Y\fn,n  and  ( Zmnt„  -  Z i„)  and 
then  the  y"n  are  set  to  the  fractional  angle  subtended  by  the  point  at  the  intersection  of 
the  centreplane  and  the  deck. 


y'mn  =  -arctan 


2 

i 

7T 


(1n 


(ZMn,n  Zmn) 


(6.4) 
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(6.5) 


where  7„  is  the  normalization  factor  for  the  nth  offset  station,  given  by 

_  Z\fn,n  ~  n 

ln  ~  n,„,n 

By  normalizing  the  Y  and  Z  values  before  calculating  the  angles,  this  method  avoids  the 
clustering  of  the  j/"-values  obtained  using  the  method  2.  Figure  17  shows  the  thin  section 
of  Figure  16  with  y"-values  as  calculated  using  the  method  3;  the  clustering  of  y"-values 
has  been  avoided.  However,  the  normalization  can  bring  back  non-local  effects;  use  of 
this  method  on  the  hull  of  Figure  14  generates  similar  kinks  in  the  y'Mines. 


4. 


A  composite  of  methods  2  and  3  can  be  defined  as  follows.  The  y'^n  are  calculated  using 
method  3  up  to  midships.  Beyond  midships,  the  normalization  factor  retains  a  constant 
value:  whatever  its  value  was  at  the  midship  station.  Thus,  if  Xj  is  the  midship  station, 
then 


7n  = 


^Mn,n  Z\n 

„  yMn,n 

zM,,i  -  Zij 

Yu,j 


for  n  <  j 
for  n  >  j 


(6.6) 


This  method  works  well  for  most  destroyer  and  frigate  hulls  in  which  the  problems  in  the 
method  2  tend  to  occur  near  the  bow,  while  the  problems  in  method  3  tend  to  occur  near 
the  stern.  As  with  methods  2  and  3,  it  will  not  work  with  offsets  with  negative  Y-values. 


5.  The  y"  n  are  calculated  in  a  user-supplied  sub-program.  For  complex  hulls,  it  is  possible 
that  none  of  the  above  methods  will  work  satisfactorily.  In  this  case  the  user  may  deter¬ 
mine  the  y"- values  by  some  other  method  of  his  choosing,  but  must  supply  a  sub-program 
to  perform  the  calculation. 


6.  The  y"  n  are  read  from  a  user-supplied  file. 


As  has  been  illustrated  above,  the  first  four  methods  for  determining  y"  offered  by  HLLSRF 
all  have  disadvantages  for  certain  hull  forms.  Other  methods  for  defining  y"  have  been  proposed, 
notably,  conformal  mapping  onto  a  unit  circle22,23  and  requiring  that  the  coordinate  system 
be  orthogonal24.  The  conformal  mapping  method  suffers  from  the  same  non-local  behaviour 
as  methods  1  and  3.  The  method  of  orthogonal  coordinates  is  not  always  well-defined  if  there 
are  hull  discontinuities  and  will  not  always  meet  the  requirement  that  the  boundaries  of  the 
spline  surface  S  be  curves  of  constant  y".  Thus,  to  date,  no  method  for  defining  y'^n  has  been 
discovered  that  is  fully  satisfactory  for  all  typical  hull  forms.  Therefore,  HLLSRF  offers  the 
user  the  two  additional  options,  5  and  6,  whereby  the  method  for  determining  y"n  can  be 
tailored  to  the  particular  hull  for  which  a  representation  is  desired.  As  the  complexity  of  the 
given  hull  increases,  so  also  increases  the  likelihood  that  the  user  will  have  to  resort  to  one  of 
these  options.  Development  of  more  satisfactory  methods  for  generating  the  y"n  is  an  area 
deserving  of  further  research. 
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6.2.2  Adjustment  of  the  t/'-values  at  Fully  Discontinuous  Stations 


At  a  fully  discontinuous  station,  the  offset  file  must  contain  two  sets  of  offset  points:  one 
set  which  lies  on  the  surface  S  just  in  front  of  the  discontinuity,  and  another  set  which  lies  on 
the  surface  just  beyond  the  discontinuity.  In  the  offset  file  this  is  implemented  by  having  two 
offset  stations,  x„d  and  xnd+i ,  at  the  fully  discontinuous  station:  xn<1  =  xn<j+j.  The  y"- values 
of  these  offset  points  are  y"„d  and  y"nd+ 1- 

The  discontinuity  will  extend  only  over  a  portion  of  the  station.  On  each  of  the  two  offset 
stations,  the  two  offset  points  lying  at  each  end  of  the  discontinuous  portion  of  the  station  must 
be  marked:  on  the  station  nj,  these  are  the  points  mi  and  m2;  on  the  station  nj  +  1  they  are 
the  points  m3  and  m^.  Hence  on  station  nj,  the  discontinuous  portion  is  y"  €  (j/m1nd'ym2nd)> 
while  on  the  station  n<i  +  1  the  discontinuous  portion  is  y"  €  (j/^3Rd, y"4tJd).  Outside  the  region 
of  discontinuity,  the  two  stations  are  identical:  hence,  the  offset  points  are  also  required  to  be 
identical  outside  the  region  of  discontinuity.  Thus, 


mi  = 

m3 

(6.7) 

(Fm„d, 

Zmnd )  = 

{Ym,nd+\,Zm,nd. fl)  for  m  <  TTlj 

(6.8) 

(Ym„d, 

Zmnd )  = 

(lm+mi  —  m2  ,nd+l  i  Zm+m4  _ m2  ,nd+l )  for  m  >  771 2 

(6.9) 

Figure  18  shows  the  offset  points  on  the  two  stations  used  to  define  the  end  of  the  propeller 
bossing  in  Figure  9.  The  points  on  station  xn  are  marked  with  dots,  while  those  on  station  xn+i 
are  marked  with  crosses.  The  centreplane  is  also  shown.  In  this  case  the  region  of  discontinuity 
extends  all  the  way  to  the  keel  so  that  mi  =  m3  =  0. 

So  that  the  hull  splines  are  continuous  outside  the  range  of  discontinuity,  it  is  also  necessary 
that  the  y-values  of  the  offset  points  on  the  two  stations  are  identical  outside  the  region  of 
discontinuity.  However,  when  the  y"-values  are  calculated  using  methods  2  or  3,  this  will 
usually  not  be  the  case.  For  method  1,  the  mismatch  is  caused  by  the  the  normalization  by  the 
total  arclength,  since  the  two  stations  will  usually  have  different  total  arclength.  For  method  3, 
it  is  caused  by  the  fact  that  the  Y-values  at  the  deck  edge  may  differ  for  the  twro  stations,  as 
may  the  Z-values  at  the  keel;  thus  the  normalization  of  the  angles  causes  differing  y"-values  for 
the  two  stations. 


Hence,  it  is  necessary  to  adjust  these  values  to  y'mu  for  which 


y'mnd  =  y'm,nd+l  for  m  <  ml  '  (6-10) 

y'mnd  =  J/m+m«-m2,nd+ 1  for  m  >  m2  (6.11) 


This  is  done  by  defining  a  transformation  g(y")  which  is  applied  to  the  second  station  of  offsets 
and  all  subsequent  stations. 


with  the  condition  that 


9(y 

g(y 


Vmn 

Vmn 

for  n  <  nj 

(6.12) 

Vmn 

9{ymn 

)  for  n  >  nd 

(6.13) 

"  \ 

~  ^m,nd 

for  m  <  m\ 

(6.14) 

"  \ 
m,nd+l  > 

=  y"  m 

+m2tnd  for  m  >  m2 

(6.15) 
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and  {t^}  are  constructed  with  multiple  knots  reflecting  the  locations  of  the  knuckles 
and  discontinuous  stations.  For  example,  there  is  a  knot  of  multiplicity  kv  -  1  in  the  j/-knot 
sequence  at  each  of  the  knuckle  y- values,  st. 

In  order  to  determine  the  hull  spline  coefficients  ajn  and  /?jn,  equations  (3.3)  and  (3.4) 
are  first  rewritten  as  follows. 


Nv 

Y  =  B^N(x)Bhkv(y) 

(6.21) 

j= i 

N„ 

Z  =  D^Vj{x)Bi>kv(y) 

(6.22) 

i= i 

Nx 

Vj(x)  =  Yi  QjnBn,kx(x) 

(6.23) 

n=  1 

Nx 

uj(x)  =  0jnBn,kx(x) 

n=l 

(6.24) 

The  specification  of  the  two-dimensional  spline  coefficients,  ajn  and  (3jn,  may  then  proceed  by 
fitting  a  series  of  data  sets  with  one-dimensional  splines. 

First,  at  each  offset  station,  the  data  set  {(ymn,Fmn)  :  m  =  1, . . . ,  A/„}  is  fitted  according 
to  equation  (6.21)  to  yield  /ij(x„).  Similarly,  the  data  set  {(ymn,Zmn)  :  m  =  1  is 

fitted  according  to  equation  (6.22)  to  yield  i/j(xn).  Then,  for  each  j  =  1,...,JV„,  the  data 
set  {(xn,/Zj(xn))  :  n  =  1  is  fitted  according  to  equation  (6.23)  to  yield  aJtl,  and 

{( xn,t/j(xn ))  :  n  =  1,.. .  ,Na}  is  fitted  according  to  equation  (6.24)  to  yield  /?jn. 

Each  of  these  splines  is  determined  using  a  single  spline  fitting  sub-program  based  on  the 
subroutine  BSMTH  described  by  Hally26.  The  algorithm  used  has  the  following  properties. 

1.  The  spline  may  be  smoothed  by  an  amount  specified  by  the  user.  ‘Errors’  may  also  be 
assigned  to  the  data  points  to  allow  the  splines  to  pass  more  closely  to  or  further  from  a 
given  point.  The  errors  may  be  changed  interactively  by  the  user. 

2.  The  user  may  assign  ‘stiffness  weights’  to  portions  of  the  spline  to  allow  or  inhibit  cur¬ 
vature  of  the  spline  in  a  given  region;  thus  a  high  stiffness  weight  is  assigned  to  the 
flat  bottom  of  a  hull  but  a  low  stiffness  weight  near  the  bilge.  The  stiffness  weights  are 
calculated  automatically  but  may  be  altered  interactively  by  the  user. 

3.  The  spline  knots  are  independent  of  the  offset  points;  moreover,  unlike  most  spline  al¬ 
gorithms,  the  number  of  B-splines  comprising  the  spline  can  exceed  the  number  of  data 
points.  This  allows  flexibility  in  the  choice  of  knots  and  great  freedom  when  digitizing 
offset  diagrams  to  create  the  offset  data  file. 

4.  After  each  spline  has  been  calculated,  it  may  be  plotted  and  altered  interactively  by  the 
user  by  varying  the  required  amount  of  smoothing,  the  data  point  errors,  and  the  stiffness 
weights. 


33 


6.4  Calculating  the  Subsidiary  Splines 


Once  the  hull  splines  have  been  calculated,  the  subsidiary  splines  for  the  deck  edge,  the  wa¬ 
terline,  the  bow  profile,  the  stern  profile,  and  the  curve  of  sectional  areas  are  determined.  Each 
of  these  is  optional.  The  subsidiary  spline  curves  are  all  calculated  using  the  same  interactive 
spline  sub-program  (based  on  BSMTH)  that  is  used  for  the  hull  splines. 

6.4.1  Calculating  the  Deck  Edge  Splines 

The  OFFSRF  representation  of  offset  data17  allows  the  user  to  mark  those  offset  points 
which  lie  on  the  deck  edge.  On  the  station  x„,  denote  the  number  of  the  point  on  the  deck 
edge  by  mj.  Then  the  data  set  {(x„,ymd)n)  :  n  =  1,...,JV4}  can  be  splined  to  yield  th  ’eck 
edge  curve  yj(x).  The  knot  sequence  for  this  spline  is  constructed  so  that  the  spline  allows 
discontinuity  of  first  derivative  at  all  deck  edge  corners,  and  discontinuity  of  the  spline  itself  at 
steps. 


6.4.2  Calculating  the  Waterline  Spline 

The  user  may  set  the  design  waterline,  Zjw,  to  any  desired  height  above  the  baseline 
provided  Zdw  <  D.  At  a  series  of  stations  {x^},  the  equation 

fz{xi,Vi)  =  Zdw  (6.25) 

is  solved  for  y,  using  the  Newton- Raphson  method.  The  resulting  data  set  {(x;,y,)}  is  then 
splined  to  yield  the  waterline  curve  yw{x). 

6.4.3  Calculating  the  Bow  and  Stern  Profile  Splines 

The  bow  profile  spline  is  calculated  by  starting  on  the  keel  and  searching  for  a  series  of 
points  for  which 

fv(xi,yi)  =  0  (6.26) 

The  search  is  started  by  choosing  a  series  of  values  for  x,  and  solving  for  y<.  If  the  solution 
procedure  fails  (as  will  happen,  for  example,  at  the  forward  station  of  a  bulbous  bow)  then  a 
series  of  values  i u  is  chosen  and  equation  (6.26)  is  solved  for  x,.  This  procedure  of  alternating 
between  solving  for  x  and  solving  for  y  is  continued  until  the  deck  edge  is  reached.  The 
resulting  series  of  points  {(x,  ,  y,)}  is  then  split  into  sub-series  over  which  x,  varies  monotonically: 
for  example,  {(1.0, 0.1), (2. 0,0. 2), (3.0, 0.2), (1.0, 0.2)}  would  be  split  into  the  two  sub-series 
{(1.0, 0.1), (2.0, 0.2), (3.0, 0.2)}  and  {(3. 0,0. 2), (1.0, 0.2)}.  The  sub-series  are  then  data  sets 
which  can  be  splined  (perhaps  after  re-ordering  to  ensure  that  x,  is  increasing)  to  yield  the 
different  segments  of  the  bow  profile  curve  y™(x).  The  procedure  to  calculate  the  stern  profile 
is  similar. 
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7  Concluding  Remarks 


The  designer  of  any  useful  hull  representation  system  is  beset  by  two  main  obstacles:  the 
complexity  and  variety  of  hull  forms,  and  the  efficiency  and  ease  of  use  demanded  by  the  user. 
HLLSRF  has  been  designed  to  surmount  both  obstacles  in  a  manner  which  will  find  favour 
with  hydrodynamicists. 

Efficiency  and  ease  of  use  has  been  addressed  by  using  the  station  number  as  the  parameter 
x :  ease  of  use  is  provided  as  this  is  the  parameter  of  choice  for  nearly  all  hydrodynamic 
computations;  efficiency  is  provided  by  avoiding  costly  spline  inversions  in  favour  of  efficient 
spline  evaluations.  These  objectives  have  also  been  addressed  in  the  implementation  of  the 
computer  code  itself: 

•  data  files  are  structured  so  that  they  are  easy  to  read  and  easy  to  edit, 

•  sub-programs  are  provided  in  the  HLLSRF  library  which  perform  many  commonly  used 

calculations,  and 

•  the  code  is  modular,  efficient,  well-documented,  and  provides  clear  error-reporting. 

The  ability  to  describe  complex  hull  shapes  has  been  addressed  by  allowing  the  parametric 
domain,  £>,  to  be  non-rectangular.  Figures  2,  6,  9,  8,  and  15  are  an  indication  of  the  variety  of 
hull  shapes  that  HLLSRF  can  represent. 

The  program  HLLSPL  has  been  provided  for  generating  HLLSRF  representations  from 
offset  data.  It  must  be  recognized  that,  while  every  effort  has  been  made  to  make  the  program 
HLLSPL  easy  to  use,  the  level  of  effort  required  to  generate  a  hull  representation  will  always 
rise  as  its  complexity  increases.  Easing  the  generation  of  complex  hulls  is  the  area  in  which 
future  research  effort  can  be  most  profitably  spent. 

HLLSRF  differs  from  most  hull  representation  schemes  in  that  it  was  developed  to  provide 
a  method  of  representing  existing  hulls  accurately  so  that  their  hydrodynamic  properties  can  be 
calculated.  Similar  representation  schemes  have  been  developed  for  the  purposes  of  designing 
new  hulls8-16  but  Hally14  has  pointed  out  that  the  requirements  for  these  two  objectives  can  be 
quite  different.  To  date,  no  representation  scheme  has  been  designed  that  is  fully  satisfactory 
for  both  purposes.  When  a  user  wishes  to  represent  a  hull  for  which  inadequate  offset  data  have 
been  provided,  a  problem  arises:  the  hull  generation  task  is  one  partly  of  approximation  to  the 
offsets  and  partly  of  design  of  those  portions  of  the  hull  which  are  not  described  adequately 
by  the  offsets.  All  current  systems,  including  HLLSRF,  require  some  user  ‘artistry’  when 
applied  to  such  hulls.  However,  it  is  also  important  to  note  that,  while  HLLSRF  is  not  yet  the 
perfect  hull  representation  system,  it  is  currently  at  a  state  of  development  which  surpasses 
the  programs  in  which  it  is  to  be  applied.  That  is,  if  a  hull  is  sufficiently  complex  that  it  is 
difficult  to  represent  properly,  then  it  is  probably  too  complex  to  be  analyzed  satisfactorily 
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